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Abstract. We present an analysis of one-dimensional models of dynamical 
systems that possess "coherent structures"; global structures that disperse 
more slowly than local trajectory separation. We study cocycles generated 
by expanding interval maps and the rates of decay for functions of bounded 
variation under the action of the associated Perron-Frobenius cocycles. 

We prove that when the generators are piecewise affine and share a common 
Markov partition, the Lyapunov spectrum of the Perron-Frobenius cocycle has 
at most finitely many isolated points. Moreover, we develop a strengthened 
version of the Multiplicative Ergodic Theorem for non-invertible matrices and 
construct an invariant splitting into Oseledets subspaces. 

We detail examples of cocycles of expanding maps with isolated Lyapunov 
spectrum and calculate the Oseledets subspaces, which lead to an identification 
of the underlying coherent structures. 

Our constructions generalise the notions of almost-invariant and almost- 
cyclic sets to non-autonomous dynamical systems and provide a new ensemble- 
based formalism for coherent structures in one-dimensional non-autonomous 
dynamics. 



1. Introduction 

Transport and mixing processes play an important role in many natural phenom- 
ena and their mathematical analysis has received considerable attention in the last 
two decades. The geometric approach to transport includes the study of invariant 
manifolds, which may act as barriers to particle transport and inhibit mixing. So- 
called Lagrangian coherent structures were introduced ( [HYOOllHOl] ) as finite-time 
proxies for invariant manifolds in non-autonomous settings. The ergodic-theoretic 
approach to transport includes the study of relaxation of initial ensemble densities 
to an invariant density, with a special focus on initial densities that relax more 
slowly than suggested by the rate of local trajectory separation. Such slowly decay- 
ing ensembles have been studied as "strange eigenmodes" ( [LH041 IPP031 IPPE07] in 
fluids and have been used to identify almost- invariant sets |D J99|. lF05l IFP08| [F08] ) . 
Until now, a suitable framework for the ergodic-theoretic approach that deals with 
truly non-autonomous dynamics has been lacking. The main aim of this work is to 
develop the fundamental structures and results that will support a non-autonomous 
theory for an ensemble-based approach to coherent structures. 

We study non-autonomous one-dimensional dynamical systems that are given by 
compositions of expanding interval maps, and their action on ensembles represented 
by probability densities. The time evolution of a density is given by the Perron- 
Frobenius operator. For a single piecewise-expanding map these densities evolve 
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toward an equilibrium distribution which is absolutely continuous (see |LY73| ) : 
when the map is transitive, this equilibrium distribution is also unique. Thus the 
equilibrium distribution is an eigenfunction of the Perron-Frobenius operator with 
eigenvalue 1. The exponential rate of convergence to equilibrium is governed by the 
spectrum of the Perron-Frobenius operator. When restricted to the space of func- 
tions of bounded variation (BV), the Perron-Frobenius operator is quasicompact 
(see |HK82j ). meaning that there are only finitely many spectral points of modu- 
lus greater than the essential spectral radius, and each is an isolated eigenvalue of 
finite multiplicity. We will say an eigenvalue is exceptional if it is different from 
1 and has modulus greater than the essential spectral radius. It is known that in 
the BV setting the essential spectral radius is determined by the long-term rate 
of separation of nearby trajectories. Eigcnfunctions corresponding to exceptional 
eigenvalues relax more slowly to equilibrium than suggested by the local separation 
of trajectories. The existence of such eigenf unctions has been attributed to the 
presence of "almost- invariant sets" (see jDJ99[ iDFSOOl IF07] ). 

Exceptional eigenvalues have previously been found by considering piecewise- 
affine expanding maps with a Markov partition ( B96 *, DFSOO], KR04]). When 
restricted to the space of step-functions constant on the Markov partition intervals, 
the associated Perron-Frobenius operator becomes a finite dimensional operator. 
In the present work we extend these results to the non-autonomous setting. Instead 
of iterating a single map, we consider a cocycle of maps and its associated Perron- 
Frobenius cocycle. The appropriate way to describe exponential rate of convergence 
to equilibrium is via the Lyapunov spectrum of the Perron-Frobenius cocycle. As 
the Perron-Frobenius operator is a Markov operator, the Lyapunov spectrum is 
contained in the interval [— oo,0]. We look for exceptional Lyapunov exponents, 
namely those greater than the essential upper bound of the Lyapunov spectrum 
but less than zero. 

We obtain a Lyapunov spectral decomposition for the Perron-Frobenius cocycle 
into invariant subspaces with given Lyapunov exponents (see Corollary |4|) . This 
relies on a new version of the Multiplicative Ergodic Theorem (see Theorem , 
which provides an invariant splitting into Oseledets spaces even when the genera- 
tors are non-invertible. Our new version strengthens the standard Multiplicative 
Ergodic Theorem (see, for example [A98, Theorem 3.4.1]) where only an invariant 
flag of nested subspaces is supplied. 

We demonstrate the existence of slow-mixing coherent structures by constructing 
periodic (see Theorem [TJ and non-periodic (see Theorem [5]) examples of Lebesgue 
measure-preserving one-sided cocycles with exceptional Lyapunov exponents. In 
each case, we calculate algebraically the Oseledets subspaces associated with the 
largest exceptional exponent and verify that the second largest Oseledets space 
captures the coherent structures. 

Finally, we present an algorithm for approximating the Oseledets splitting, which 
is based on a new computational approach suggested by the proof of Theorem 3. We 
demonstrate the effectiveness of the algorithm, by approximating some Oseledets 
subspaces numerically. 
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2. Preliminaries 

We study the Perron-Frobenius operator of compositions of expanding maps. We 
first introduce the necessary notation and relevant results for autonomous systems, 
and then extend this to the non-autonomous case. 

2.1. Autonomous systems. We say that T : I — ► I, where / = [0, 1] or S , is an 
expanding map if there exists a finite partition oq = < a% < . . . < ajj = 1 such 
that, for each i = 1, . . . , M, T is continuous on (aj_i, aA and extends to a C 2 map 
on [flk_i,Oi] satisfying |DT|( 04 _ 1)0i) | > 1. 

The Perron-Frobenius operator for an expanding map T : I — ► I is defined, for 
an L 1 function / : I — > R, by 

y£T~ 1 (x) 

In LY73 , the Perron-Frobenius operator is used to prove that expanding maps 
have an absolutely continuous invariant probability measure. The key step of their 
proof is to show that the Perron-Frobenius operator contracts the norm on a suit- 
able space of functions: the functions of bounded variation. 

The variation of a function / : / — > R on a subinterval A C / is defined by 

var/ := var f{x) = sup < \f(xi) - f(xi-i)\ : Xi £ A, x < ■ ■ ■ < x k \ . 

Given / G L°° C L 1 , the variation is defined by var/ / = infjvar/ g ; f = g a.e. }. 
We denote by BV the Banach space 



BV = {/ S i°° : var / < ooj 



equipped with the norm ||/|| = max{||/||ii, var/ /}. We denote Lebesgue measure 
on I by m, and / € BV is called a (probability) density if / > on / (and 
H/Hii = 1). We write BV+ = {/ e BV : f > 0}. 

The Perron-Frobenius operator is Markov: that is, if / € i 1 is a density, then 
V '/ is also a density and H'P/Hl 1 — II/IIl 1 - A probability density /* satisfying 
V f* = f* is an invariant probability density for T. 

Keller |K84j shows that the Perron-Frobenius operator of an expanding map has 
at most finitely many exceptional eigenvalues. 

Theorem (Keller (1984)). Given an expanding interval map T : I —> I, its 

Perron-Frobenius operator V acting on BV has essential spectral radius 

l/n 



lim sup 



n — >oc , 



1 



> x £\\D(T»)(x)\, 

and all other spectral points are isolated eigenvalues of finite multiplicity. 

Exceptional eigenvalues have a distinguished dynamical significance as their 
eigenf unctions are associated with relaxation to equilibrium at exponential rates 
slower than the rate suggested by the average local separation of trajectories 9. For 
example, if Vg — Xg with 6 < |A| < 1 then an initial density /* + ag, a / will 
relax to /* at a rate slower than 9. 

Dellnitz and Junge [DJ99 suggested that positive real Perron-Frobenius eigen- 
values near to 1 correspond to almost-invariant sets] more precisely, they suggested 
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the sets A + := {g > 0} and A~ := {g < 0} formed an almost-invariant partition 
of the state space. Dellnitz et al. |DFS00j showed the converse, presenting a class 
of interval maps with almost-invariant sets and proving the existence of excep- 
tional eigenvalues. Froyland |F07j constructed a two-dimensional hyperbolic map 
with almost-invariant sets and proved the existence of an exceptional eigenvalue. 
Numerical methods have been developed ( |DJ99( F05. F08]) for the computation 
of exceptional eigenfunctions and almost-invariant sets; these have been applied 
successfully in molecular dynamics QSHD99 ), astrodynamics ( |D+05| ), and ocean 
circulation f |F+07| ). 

Our intent in the present work is to generalise the notion of almost-invariant sets 
in autonomous systems to that of coherent structures in non-autonomous systems. 
The latter will represent structures that are perhaps quite mobile, but disperse at 
rates slower than suggested by local trajectory separation. 



2.2. Non-autonomous systems. We will examine exceptional spectral points in 
the non-autonomous case, and study compositions of expanding maps taken from 
a finite collection, and composed in order according to given sequences. 

Let a be an ergodic automorphism of a probability space (fl, 7i, p) that preserves 
the probability p. Given a measurable/topological/ vector space X, a (one-sided) 
cocycle over a is a function H : Z + x £1 x X — ► X with the properties that for all 
x iE X and lo £ fl: 

• -ff (0, lo, x) = x; 

• for all m, n £ Z + , H(m + n, u>, x) = H(m, o~ n u>, H(n, lo, x)). 

We sometimes write H( n \ui)(x) for H(n,u,x), and H(uS){x) for H{\,lo,x). The 
generator of a cocycle H is the mapping H : Q — > End(X) given by H(u>) = H(u>). 
Since the cocycle is uniquely determined by H, we occasionally refer to H itself as 
the cocycle when no confusion can occur. 

In the sequel, a will frequently be a (left) shift, defined by (au>)i = cjj+i, acting 
on a two-sided sequence space (Q,,!F,p) on K symbols, where f2 C (Zk) z , Zk = 
{1, ... ,K}, is invariant under a. The shift a preserves the probability p and is 
ergodic with respect to p. 

Definition. Let {Ti}i & z K , be a collection of expanding maps of /, and let Vi : 
BV — * BV be the Perron-Frobenius operator associated to : J — ► I. The map 
cocycle generated by {Ti}i^z K , denoted by $ : Z + x f2 x / — > I, is defined to be 
the one-sided cocycle with generator $(w) = T Wo £ {Ti}i & z K - Associated to $ is 
the Perron-Frobenius cocycle V : Z + xOx BV — > BV, which is defined to be the 
one-sided cocycle with generator V(w) — V^ Q £ {Vi}iez K - 

Notice that even though we use a two-sided shift space, we only form one-sided 
cocycles not two-sided cocycles. This is because the expanding maps are non- 
invertible, as are their Perron-Frobenius operators. 

We say a cocycle is periodic if the underlying shift space f2 is generated by a 
single element: that is, there exists R £ N, called the period, and w £ such 
that f2 = {lo, a~u>, . . . cr fl ~ 1 w}; we say a cocycle is autonomous if f2 contains a single 
element. 
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3. QUASICOMPACTNESS OF THE TRANSFER COCYCLE 

Information about the exponential decay rates of the Perron-Frobcnius cocycle 
is given by its Lyapunov spectrum. 

Definition. We denote by A(o>, /) the Lyapunov exponent of / € BV, defined 
A(w,/) = limsup-log||^")(a;)/||. 

n^oo ri 

We define the Lyapunov spectrum A(V(ui)) C M of the Perron-Frobcnius cocycle 
at lo to be the set 

A(P(«)):={A(w,/):/6BV}. 

The exponential rate of decay that can be expected purely from the local expan- 
sion is the essential upper bound 

d{io) := - inf limsupilog|L>$ ( " ) (w)(a;)| < 0, 

n — >oo ri 

which we denote by $ if independent of u>. Points in A("P(ct>)) that are greater than 
indicate the presence of large-scale structures that reduce the rate of mixing 
of the system, except for the maximal Lyapunov exponent, 0, which is associated 
with an invariant density. We refer to Lyapunov spectral points in the interval 
($(w),0) as exceptional. 

In order to find systems with exceptional Lyapunov spectrum, we restrict our 
attention to map cocycles generated by piecewise-affine maps with a Markov parti- 
tion. We say an expanding map T : I — > I is piecewise-affine if there is a partition 
a = < ai < . . . < a m = 1 such that T has constant derivative on each in- 
terval (aj_i,a.j). Recall that, for a map T : I -> I, a partition B = {B k }^f =1 of 
/ into intervals is called a Markov partition if for each pair 1 < i,j < M such 
that Bi n T(Bj) ^ 0, we have Bi C T(Bj). Associated to T is a transition matrix 
T = (7i,j)i<i,j<M, where jij = 1 if T(Bj) D Bi and otherwise. Given a partition, 
we denote the set of partitioning points by Sjg = I\ Use/3 B. 

For a Markov partition B, we let x(B) denote the space of step-functions L — > M 
that are constant on the intervals of B. We say B is a common Markov partition 
for a collection of maps {Ti}nzz K if it is a Markov partition for each map Tj. If the 
generators of the map cocycle have a common Markov partition B, then x(B) C BV 
is an invariant subspace for the associated Perron-Frobenius cocycle. 

We now fix a collection T := {T{\i e z K of piecewise-affine expanding maps 
with a common Markov partition B. Clearly y(£>) is an invariant subspace for 
the Perron-Frobenius cocycle. Let F be the quotient space BV/x(£>) with norm 
WIWf ■= ||/ - Qf\\, where Q : BV — > x(#) is the projection 

Qf(x) = <*B(f)XB(x), a B (f) = ± [ f(s) ds. 

We identify an element / + x(B) <E F with the function / — Qf. Thus / S F is 
characterised as having zero mean on each interval B & B: that is, for each B e B, 
J B f(s)ds = 0. 

The following Lemma bounds the growth rate for functions in F; we will shortly 
see that ${oj) is the essential upper bound for the Lyapunov spectrum of V(u>). 
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Lemma 1. Let $ be a map cocycle generated by piecewise-affine expanding maps, 
with a common Markov partition B. Then \\V^(uj)\f\\ < sup 7 ^ D ^)^ f or every 
n £ N, and thus 

supA(w,/)<0(w). 
feF 

Proof. The common Markov partition B allows for the following estimate of the 
variation of the Perron-Frobenius cocycle: 

var7>< n >(w)/ < ^ var( XA • V in) (cu)f) 
AeB 

f(y) \ 



53 var X a(x) ■ ^ 

AeB V «£[*<»> Ml- 1 (a; 



ye[*(")( W )]- 1 (x) 1 V 

v / v XA(^ (n) H(y)) •/(!/) ' 

AGS \j/e[*(«)(w)]- 1 (x) 1 v 
AGS V b 6 BW( U ) ' V U f 



AeB see<"»(") 

* ( "'(w)B=A 



XB 



|£>*( n )(w)| 



where = V/™= $ w (w) _1 £- Thus 

Xb-/ 



(2) y&rV {n \w)f < ^ var ( 



|D*(»)(w) 



We show that var/(/ — Q/) > || / — Q/Hl 1 for / £ BV, from which it follows that 
= var/(/ — Qf). 

\\f-Qf\W = Y. I \(f-Qf)(s)\ds 

BeB jB 

< V \B\ (esssup(/ - Qf)(x) - essinf(/ - Qf)(x)) 
BeB V xeB xes y 

< ^ fcsssup(/ - Qf)(x) - essinf(/ - Q/)(ar)) 

< 53 var(/ - Qf) 



B 

BeB 

(3) <var(/-Q/). 
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Let t;(f(x)) = | lim y y x f{y) — lim^x^ f(y)\ denote the jump of / at x. Now, 
suppose / G F. Applying ([2]) we have 

XB-f 



\\v {n) (")f\\ < E v f 



D$(«) M| 



\D<$> 

/ 



< 3var 



/ |£»$(™)(cj)|' 
. 3 . 

(4) = T i^")Mr ll/ "' 

since / G F by (J3j> , giving the first part. Hence, for any / G F, 

and the second result follows. □ 

We now prove that the exceptional Lyapunov spectrum of V is contained in 
A("P(w)| x (g)). For the autonomous case, see for example |BK981 Lemma 3.1]. 

Proposition. Let $ be a map cocycle generated by piecewise-affine expanding maps, 
with a common Markov partition B. Then the Lyapunov spectrum of the Perron- 
Frobenius cocycle V(u>) satisfies 

A(P( W ))C{.<%)}UA(P( W )| #) ), 

and thus V(uS) has at most #B exceptional Lyapunov exponents. 

Proof. Each / G BV has a unique decomposition as a sum / = f x rm + /f, where 
f x (B) ■= Q(f) G X(B) and f F := / - Q(f) <E F = BV/x(B). Notice that 
A(w, /) < max{A(w, / x (b)), A(w, / f )}. Thus, if X(w, f x (B)) < A(w, /f), then we have 
A(w,/) < by Lemma U Otherwise, we have X(u>,f x ^)) > A(w,/ F ), in which 
case A(w, /) = X(uj,f x{B) ). □ 

4. A STRONGER MULTIPLICATIVE ERGODIC THEOREM FOR NON-INVERTIBLE 

MATRICES 

By Proposition [3j for each uj G 17, all exceptional Lyapunov exponents of V(lu) 
are contained in the Lyapunov spectrum A(V (oj)\ x (b)) ■ We now represent V(uj)\ x (b) 
as a matrix cocycle. 

The set {x(Bi)}iLi forms a basis for and thus each / G x(^0 ma Y be 

written as / = 'Y2^L\ v iXB i in a unique way. Similarly, given v G K M , we write 
(u) := Yli=i v iXBi for the corresponding function in BV. 

For T G T, the matrix P = (pij)i<i.j<M, where 

^ = Ml-(^ 1< U <M, 
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represents the Perron-Frobenius operator for T with respect to the basis {x(Bi)}f£ 1 
of x(&) ( see ) f° r example, }BG97|, p. 176]). That is, for each v £ R M we have 

V (v) = (Pv) . 

Let Pi denote the matrix representing the restricted Perron Frobenius operator 
Vi\ x (p) with respect to the basis {x{Bi)}iLv The matrix cocycle A : Z + xf2xR M — > 
R M is the one-sided cocycle with generator A(oj) — P Wo . 

Thus for each oj G O, all exceptional Lyapunov exponents of V(ui) are captured 
by the Lyapunov spectrum of the cocycle A(A(uj)) = A(V '( w )lx(B)) ■ 

The Multiplicative Ergodic Theorem for one-sided cocycles (see, for example, 
IA981 Theorem 3.4.1]) provides us with a description of the asymptotic behaviour 
of the matrix cocycle A(oj). It reveals that the Lyapunov spectra = Ai > A2 > 
■ ■ ■ \g > —00 of A(lo) coincide for all w in a er-invariant il C £1 of full p-measure. 
Moreover, it states that for each w £ fi, a Lyapunov exponent A(w, v) of v £ 
R M for A(u>) is determined by the position of v within a flag of nested subspaces 
{0} = Vt(w) C • ■ ■ C V 2 (pj) c Vi(uj) = x(B). Specifically, for each i = l,...,£, 

(5) \{w,v) = * i <=^veV i (w)\Vi+i(w). 

In addition, the flag of subspaces is preserved by the action of the cocycle: for 
i= I,...,?, 

A(u)Vi(u) C Vi(au). 

For two-sided matrix cocycles (see, for example, |A98| Theorem 3.4.11]), by inter- 
secting the corresponding subspaces of the flags for the cocycle and for its inverse, 
one obtains an Oseledets splitting: that is, for each u) £ Q we have a decomposition 
K M = 0- =1 Wi{u>) such that for i = 1, . . . ,£, 

o6Wj(w)\{0}=^A(w,ti)=^ 

and 

A(cj)Wi(uj) = W t (auj). 

Our matrix cocycle A(cu) sits between these two extremes: the shift a is in- 
vertible, but the matrices {Pi}i & z K generating A(cj) are not. Because of the non- 
invertibility of the cocycle, we cannot use the standard approach described above 
to define an Oseledets splitting. The following new result relies on a push-forward 
approach to prove the existence of an Oseledets splitting even when the generators 
are non-invertible. 

Theorem 3. Let a be an ergodic invertible measure-preserving transformation of 
the space (Q,B,fi). Let A: — > A/^R) be a measurable family of matrices satisfying 

J log + \\A(uj)\\ d/i(w) < 00. 

Then there exist Xi > A 2 > • • • > \g > —00 and dimensions mi,...,m<, with 
mi + • • • + me — d, and a measurable family of subspaces Wj(ui) C R d such that for 
almost every to £ the following hold: 

(1) dim Wj (ui) = rrij ; 

(2) M. d = @ i j =1 W j {w); 

(3) A(ui)Wj(u>) C Wj(auj) (with equality if Aj > — 00 j; 
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(4) for all v £ Wj(w) \ {0}, one has 

lim (l/n)]og\\A(a n - 1 u)...A(u)v\\ = Xj. 

n — ►oo 

Proof. See Section!] □ 

Remark. It follows from part (iv) of Theorem 3 that we can determine the Lyapunov 
exponent for any vector v £ M. M \ {0} by 

i+i e+i 
A( W) «) = A i ^«60^)\ W k (uj), 

k—i k—i+1 

where we set Wg+i(uj) = {0} for all u> € Q. 

We now apply Theorem [3] to our matrix cocycle A(u>) induced by P(lu)\ x ^. 
Consider the part of the Lyapunov spectrum of A(ui) that is greater than i?(cj). Let 
r < £ satisfy A r +i < < A r . Thus, the part of A.(V (u>)\ x /m) strictly greater 

than $(u/) is precisely Ai > A2 > • • • > A r . It follows from Proposition [3] that 
the exceptional Lyapunov spectrum of V{u)) is precisely A2, • • • , X r . By defining 
Wi(uj) — {(v) : v £ Wi(oj)} for 1 < i < r, we transfer the splitting of M. M obtained 
from Theorem [3] into a splitting of x(B) to get the following result: 

Corollary. Let (f2, T \p, a) be a measurable sequence space, and $ : Z + xQx I —> I 
be the cocycle generated by piecewise-affine expanding maps, with a common Markov 
partition B. Then there exists a forward invariant full measure subset C fl, 
= Ai > • • • > A r > and mi, . . . , m r , satisfying m\ + • • • + m r < j^-B, such 

that for all u £ fl, 

(1) there exist subspaces Wi(u>) C BV, i — 0, . . . ,r, dim(Wi(uj)) = m^; 

(2) V(Lo)W l {Lu)^W l (aio); 

(3) (v) £ W z (lo) \ {0} =► \(u,v) = A,. 

Note that for periodic cocycles, ft = il. 

5. Construction of periodic cocycles with exceptional Lyapunov 

spectrum 

In this section we build a periodic map cocycle for which the Perron-Frobenius 
cocycle has an exceptional Lyapunov spectrum. 

In [DFSOOj individual maps are constructed for which the Perron-Frobenius oper- 
ator has exceptional eigenvalues. The construction uses so-called 'almost-invariant' 
sets. Given a map T : I — * / with an absolutely continuous invariant probability 
measure /i, a subset U C / is almost-invariant if 

m ~ ' 

For a map with an almost-invariant set U, the transfer of mass between U and 
I\U is low, and so we expect to find that a mean-zero function positive on U 
and negative on I\U decays to zero slowly. It is shown that for piecewise-affine 
Markov maps, one often obtains an almost-invariant set from the support of either 
the positive or negative part of the eigenfunction associated to the second largest 
eigenvalue of the Perron-Frobenius operator. 

For this first example, we construct a cocycle over a periodic shift space of period 
3 that has a cyclic coherent structure. More precisely, we take three maps, each 
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having a distinct interval from the partition J = {[0, 1/3], [1/3, 2/3], [2/3, 1]} of 
S 1 as an almost-invariant set. Post-composing these maps with the rotation by 
1/3, we form three new maps which we apply in sequence repeatedly, thus forming 
a periodic map cocycle In this way, each generator <t(ui) of the map cocycle 
moves the majority of the mass of one distinguished interval J(w) <E J into another 
interval J{auj) g J with some small dissipation. Thus a map J : Q — ► J specifies 
the location of our coherent structure. 

Theorem 1. There exist three piecewise-affine expanding maps T\, T2, T3 : S 1 — > S 1 
with a common Markov partition B that generate a map cocycle <f> : Z + x il x S 1 — > 
S 1 over the sequence space f2 C generated by a = 123 with the following prop- 
erties for i = 1,2, 3: 

(1) each map Tj preserves Lebesgue measure; 
(2; = log 1/3; 

f5) each finite dimensional restriction *Pi\x(B) of the Perron-Frobenius operator 

of Ti has no exceptional eigenvalues; 
(4) ~P(lo) has an exceptional Lyapunov spectrum that is independent of lu and 

satisfies 



A(P)n(«,o> = {l*p^)}. 

(5) the Oseledets subspaces W2(w) corresponding to the largest exceptional Lya- 
punov exponent exists for all uj g f2, and depends only on u>q. 

For periodic map cocycles, one can find Lyapunov spectral points from the eigen- 
values of the cyclic composition of Perron-Frobenius operators. 

Lemma 2. Consider a periodic map cocycle $ : Z + xQxI^flxI of period R. 
If 7] is an eigenvalue of the Perron-Frobenius operator V^ r \uj) , then 

^ e A(P( W) ). 

Proof. There exists a function / / € BV such that V^ R \uj)f = rjf. Hence for 
any k g N and < r < R, 

min {||P W (w)/ll} < H p(fcfi+r) (^)/ll < max {11^(0(^/11} 

and the result follows. □ 

Proof of Theorem]]^ Consider the partition J = { Ji, J2, J3} of S* 1 into the subin- 
tervals J; = \{i — l)/3, z/3]. Let $ : Z + xilx/^Slx/be the map cocycle with 
generator $(w) = T Wo , where the maps T = {Ti, T2, T3} are given by 



T,(i) = 3a; - J - + ^ (mod 1), z g B J ■■ 
where 

'6 7 6 1 3 

(,'=[3 6 5 8 

6 7 1 6 



1 i 

'9 



1 9 

J''*) 3 
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The graphs of T\,T2,Ts are shown in Figure [TJ note that, by construction, each 
map Ti largely maps the interval Ji into the interval Jj+i, taking indices modulo 
3: in fact, for j 6 1,2, 3, 

m(J i nT- 1 Ji +1 ) = 8 
m{Ji) ~ 9' 

Thus we have a coherent structure built around the family of intervals J : Q — ► J 
given by J(w) = J^ . 



Ti T 2 T 3 




Figure 1. Graphs of T 1} T 2 ,T 3 . 

Note also that each map Ti is piecewise-afhne expanding and there is a common 
Markov partition for T given by B — {Bi : i = 1, . . . , 9}. Notice that for each map 
T e T and interval B € B, DT\ B = 3, and so <&(u) = log 1/3 for each uj e Q. 

Moreover, for each map T £ T and interval B G B, the preimage T~ X B has 
precisely three components, each of one third of the length of B. Thus each T G T 
preserves Lebesgue measure, and hence each & n \uj), lu S fl and n € N, does also. 

As before, let Pj denote the matrix of the restriction Vi\ x ^) with respect to the 
basis {x(B)}i£z g - Here Pi — T^/3 is the one third scaling of the transition matrix 
, which is itself easily observed from the graph of T^. the (p, q)th entry of the 0-1 
matrix Ti is 1 if and only if the graph of Ti intersects the (p, q)th square of B x B. 
Each matrix Pi has a simple eigenvalue 1, and all other non-zero eigenvalues lie on 
the circle of radius 1/3: 

spec(Pi) = (1,-1/3, -1/3, 0,...,0) 

s P ec(P 2 ) = (1,1/3, 0,...,0) 

spec(Ps) = (1,-1/3, -1/6 ±i-y/3/6,0,...,0). 

Unlike in Theorem [5] in the following section, the maps used here cannot be ex- 
pressed as different rotations of a single map. 

We can find slowly decaying functions by examining the triple composition 
<&( 3 \a) = X3 o T2 o T\. The Perron-Frobenius operator V^(a), when restricted 
to the space x(#), can be represented by the matrix ^^(a) = P 3 P 2 Pi. We have 

s P ec(^ 3 V)) - (i,| ; (4±2n/TI),0 ) ...,oJ 

Since the cocycle is periodic, we find that the spectrum of A^(uj) is independent 
of to G SI. Applying Lemma[2]we have that A(P) has the two exceptional elements 
with approximate values 

A 2 w log 0.8153, A 3 w log 0.3699. 
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(w 2 ) 



-2 - 

Figure 2. The graph of (wz) € x(&) f° r Theorem 1. 



Moreover, these Lyapunov exponents are achieved by the corresponding eigenvec- 
tors of A^ 3 \uj). For lj = a, the space W2{a) is spanned by the second eigenvector 
W2 of the matrix A^(a) = P3P2P1, with approximate entries 

w 2 = (0.105, 0.193, 0.193, 0.008, -0.059, -0.059, -0.113, -0.134, -0.134), 

and the graph of (w 2 ) G x(&)i which spans Wi{oi), is shown in Figured For 
i = 1,2, W2(c*a) is spanned by (w^) ° where p : S 1 — > S 11 is the rotation 
p(x) =x + l/3 (mod 1). □ 

Evidence of the cyclic coherent structure is visible in the second eigcnfunction of 
the Perron-Frobenius operator. Note that J (a) = [0, 1/3] supports the majority of 
the mass of the positive part of {1V2) ■ Similarly, the distinguished interval J {a 1 a) — 
[(i — l)/3, i — 1,2, is picked up by (w 2 ) o p~ % . 



6. Construction of non-periodic cocycles with exceptional Lyapunov 

spectrum 

We now construct a non-periodic map cocycle with exceptional Lyapunov spec- 
trum. The map cocycle is generated by six maps, including T\ used in the previous 
example. The shift space is taken to be a subshift of finite type that has the 
Bernoulli shift on two symbols {Z% , 8, /1) as a factor. 

Let & C Zg be the subshift of finite type 

6 := {u £ {Z%f : Vfc e Z, (wk,w*+i) = iff E UJ = 1}, 



with transition matrix 



E — {Eij)i<ij<6 — 



/ 1 

1 

1 



1 

1 
\ 1 



1 \ 

1 

1 



1 

1 
10/ 



We let (7 : @ — ► be the left shift, and p the uniform measure on 8. Notice that 
({0, 1} Z , 8, (j.) is a factor of (6, a,p) via the mapping 



h{u)i = 



0, 
1. 



{1,2,3}, 
w< G {4,5,6}. 
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The six maps are obtained from T\ by rotations, and constructed so that 

(R s rnjjj n T-\l l+1 ) 8 

(6) mj-) = 9 for * = 1 ' 2 ' 3 ' 

(7) = 9 fOTi = 4 ' 5 ' 6 - 

From these maps we construct a map cocycle with a non-periodic coherent structure 
that is responsible for the slow decay. 

Theorem 2. There exists a collection S of six piecewise-affine expanding maps 
Si,...,Sb : S 1 — > S 1 with a common Markov partition B that generate a map 
cocycle $ : Z + x 6 x S 1 -> S 1 with the following properties for i = 1, . . . , 6: 

(1) each map Si preserves Lebesgue measure; 

(a; = log i/3; 

f3j the restricted Perron-Frobenius operator Pi| x (s) ^ as no exceptional eigen- 
values; 

(4) for each w G 6, A(7 3 (o;)) contains a unique exceptional exponent 

, 1 + V2 
I08-3— 

i/iere exists an Oseledets decomposition for all uj G 0, and £/ie Oseledets 
subspace W2(uj) depends only on uiq. 

Proof. Let p : S" 1 -> S" 1 be the rotation a; h-> x + 1/3 (mod 1) and let S : S 1 -> S* 1 
be the map given by 



5(a;) = 3a; - - + ^ (mod 1), x G Pj 
3 9 



j - 1 i 
9 ' 9 



i = i, 



where 5 = (3,4,3,7,0,6,1,0,6). The interval Ji = [0,1/3] is an almost-invariant 
subset of 5 1 , with m{J\ n S~ 1 Ji)/m(Ji) — 8/9. Let P5 be the matrix of Vs\ x {B) 
with respect to the basis The spectrum of Pg is 

spec(P s ) = ^l,ii^,0,...,0j . 

We define the collection of maps 5 = {Si}i e z e m terms of S and p: 



Si 


= p s 


5-2 


2 2 

= p b p 


s 3 


= S p 




= p 2 oS 


s 5 


= Sop 2 


s 6 


= p S p. 



The graphs of S\,. . . ,Se are shown in Figure [31 Note that the graph of Si is the 
same as that of T\ shown in Figure [TJ 

Let $ : Z + x 6 x S 1 — > S" 1 be the map cocycle with generator $(w) = 5 Wo G 5. 
Let J7" = {Ji}f=iJ where Ji = [(« — l)/3, i/3]. As a consequence of (j6|) and 0, we 
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Figure 3. Graphs of Si, . . . , Sq. 



have a coherent structure built around the family of intervals J : 6 — > J, where 

j{w) = [ J r 

V ^ \ Juo-3, if W > 3. 

Let Vi be the Perron-Frobenius operator of Si. Let V : Z+ x 6 x S 1 — > S 1 
the Perron-Frobenius cocycle associated to Let Pi be the matrix representing 
Vi\ x (B) with respect to the basis x(&) an d let A : Z + x Q x S 1 — > S 1 be the matrix 
cocycle with generator A(u>) = P WQ . Let R denote the matrix with Rij = 1 if 
i — j = 3 (mod 9) and otherwise. Note that R 3 is the identity matrix. For 
i = 1, ... ,6, the formula for Pi is obtained directly from the formula for Si by 
replacing p by R and replacing S by Pg. Thus, for i = 1, . . . , 6, we may write 
P = R h .P s .R r % where I = (1,2,0,2,0,1) and r = (0,2,1,0,2,1). 
One may confirm that 



spec(Pj) 



(1,-1/3,-1/3,0,. ..,0), if i < 3; 
(1,0,. ..,0), ifi>3, 



and so no map in S has exceptional eigenvalues. 

Note that whenever E t j = 1, we find k +rj =0 (mod 3). Hence for any we6, 
we have that 



A^(uj) = ^(""-l) o (P s )" o P r ( w ° 



with all inner R factors cancelling. 
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(w 2 ) 



Figure 4. The graph of (wz) G x(&) f° r Theorem 2. 

Hence for any v G K M , 

\\A^(cj)v\\ =||i^n-i .(Ps) n .R r «°v\\ 
= \\(Ps) n .R rM v\\ 
= \\(PsT-v'\\, 

where v' — R r ( u °>v. So A(A) is precisely the set of logarithms of the eigenvalues of 
Ps, and in particular, is independent of to. Thus, A('P) has a unique exceptional 
exponent log(l + \/2)/3 with approximate value log 0.8047 for every w£9. 

Let W2 be an eigenvector of Ps corresponding to the second largest eigenvalue 
(1 + \/2)/3. The graph of (w 2 ) € x(S), which 

spans W2(a), is shown in Figure 
21 Moreover, we have an Oseledets splitting for every u> E 6: for each u> £ 0, the 
function (i?~ r ( WQ )u! 2 ) spans the Oseledets subspace Wa(w) associated to log(l + 
\/2)/3 and thus W2(w) depends only onwj. □ 

As in the periodic example, the coherent structure responsible for the slow decay 
is detected by the second eigenfunction of the Perron-Frobenius operator. When 
(Jq = 1, J(u>) = [0,1/3] is the distinguished interval for $(w), and this interval 
supports the majority of the mass of the positive part of the function {w^) spanning 
W2(w). More generally, for u & Q, the positive part of {1V2) ° p~ r ^ u °^ = (R~ r ( Wo ^W2) 
is supported approximately on the interval J(u>). 

7. Numerical approximation of Oseledets subspaces 

In this section we outline a numerical algorithm to approximate the Wi(u>) sub- 
spaces. The Oseledets splittings for the cocycles in Theorem 1 and Theorem 2 
were explicitly constructed as eigenvectors. In general, the Oseledets splittings are 
difficult to compute. The algorithm is based on the push-forward limit argument 
developed in the proof of Theorem [3J 

Algorithm (Approximation of the Oseledets subspaces Wi(ui) at lu 6 fi.). 
(1) Choose M,N > and form 

^ M Ha- N u>) := (A( A V^) T ^ (M V^)) 1/2M 
as an approximation to the standard limiting matrix 

/ , . , \ 1/2M 

B(<j- n oj):= lim ( (a~ N lj) t (a~ N lj)) 
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Figure 5. Graphs of T 4 , T 5 and T e . 



appearing in the Multiplicative Ergodic Theorem. 

(2) Calculate the orthonormal eigenspace decomposition of ^^ M \a~ N uj), de- 
noted by ul M \a- N uj), i = l,...,£. 

(3) DefineWl_ M ' N \uj) := A^ N ^(a~ N ui)U < > M ^{a~ N uS) via the push forward under 
the matrix cocycle. 

(4) W^ M ' N \uj) is our numerical approximation to W%{w). 

Remarks. 

(1) Theorem 0] states that w}°°' N) (w) -> Wi(u) as N -> oo. 

(2) This algorithm also provides an efficient numerical method for calculating 
the Oseledets subspaces for two-sided linear cocycles. 

(3) There is freedom in the choice of relative sizes of M and N: in order to 
sample equal numbers of positive and negative terms of u>, we take M = 2N. 

Example. To illustrate this technique, we calculate the Oseledets subspaces W2{u>), 
oj G 0, for a non-periodic map cocycle, created from the maps of Theorem[T]and the 
sequence space 9 of Theorem [5J Unlike the example of Theorem 2, this example 
does not have a simple structure that makes it possible to relate the Oseledets 
subspaces to those of a single autonomous transformation. 

Let T = {Ti}? =1 denote the collection of piecewise-affine expanding maps of the 
circle consisting of the three maps T\,T2,T$ defined in Theorem Q] and the three 
maps T4 = poT\, T5 = p o T2 and Tg = p o T3, where p : S 1 — -> S 1 is the rotation 
p(x) = x + 1/3 (mod 1) as before. The graphs of the maps in T are shown in 
Figures Q] and [5] Let $ : Z + x O x S 1 — > S* 1 be the map cocycle generated by 
T. The collection T has a common Markov partition S = {[(«— l)/9,z/9] : i = 
1, . . . ,9}. We expect to find an exceptional Lyapunov spectrum since the cocycle 
has a coherent structure similar to that of Theorem 2, built around the family of 
intervals J : Q — * J given by 



4„, if w < 3; 
■>o;o-3) if w > 3. 



We generate a test sequence in O as follows. Let a* € {0, 1} N be the fractional 
part of the binary expansion of it: 



a* = (0,0,1,0,0,1,0,0,0,0,1,1,1,1,1,1,0,1,1,0,1,0,1,0,1,0,0,...), 
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and extend it to a two-sided sequence a * € {0, 1} Z by defining a* = for i < 0. 
We define uj* = h~ 1 (a 120 a*), where h is the 3-to-l factor defined in Section 6, and 
we take the inverse branch with u>q = 1. Note that u*6 6 has the form 

uj* = (...,1,2,3, 1,2,3, 4,... ,5, 4, 6,2,3, 1,5, 4,3, 1,5,1, 5, 4, 6,2, 6, 5,1,...), 

where the zeroth term is underlined. 

As before, we denote by Pj the matrix representation of the Perron-Frobcnius 
operator Vi\ x (B) of Tj, i — 1, ... ,6, with respect the basis x(^); an d denote by 
A : Z + x 6 x S 1 ^ S 1 the matrix cocycle with the generator A(u>) = P Wo . The 
Multiplicative Ergodic Theorem states that for almost every cj, ty( M >(cj) converges 
to a limit B(ui) as M — ► oo, and moreover 

A(A) = logspec(P). 

Calculating ^/( M )( Ci ;*) for M = 40, we find that ^( M )(w*) has a simple eigenvalue 
A2 ~ 0.81, suggesting that that V has exceptional Lyapunov exponent approxi- 
mately equal to log 0.81. 

In order to approximate the Oseledets subspace W%{u)*} numerically, we set 
M = 2N = 40, form the matrix ^ 2N \(j- N uj*) and denote by u ( ^ N) {a' N uj*) the 
eigenvector corresponding to the eigenvalue A2 . We then calculate 

A^ N \a- N uj*)u < i N \a- N ^) 

and normalize to give the vector w^ N ' N \lu*). The unit vector w^ N ' N \lu*) is our 
approximation to a unit vector spanning the subspace W^t^*)- 

Although Theorem [3] holds only for a full p- measure subset of 9, and so can tell 
us nothing about a particular sequence such as u>* , we can still check whether its 
conclusions hold in this case. Taking N — 20, we calculate for k = 0, . . . , 7, a vector 



w. 



(<r k uj*) spanning W^ N ' N) (a k w*) (see Figure© 







Figure 6. The graph of (u^'^V^*)) for fc = 0, . . . , 7. 

Recall that {ui*} J =Q = {1,5, 1, 5, 4, 6, 2, 6}. For k = 0, . . . , 7, by examining Figure 
IH and comparing with the list ( J(a k uj*))J =0 given by 

([0, 1/3], [1/3, 2/3], [0, 1/3], [1/3, 2/3], [0, 1/3], [2/3, 1], [1/3, 2/3], [2/3, 1]), 
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we see that the interval J(a k ui*) is approximately picked up by the support of the 

c (2N,N), k 

of w y 2 '(a i 
i check propert; 
is close to W^ N ' N \a k+l u)*), we calculate the quantity 



positive part of w^ N ' N '(a k u)*). 

In order to check property (iii) of Theorem 3, that is, whether A{lj*)W^ 2N ' N \uj*) 



w 



2 )± 



A(uj*)w^ N ' N \uj*) 
\{A(u*)w^ N) (w*)) \\ L1 



for N = 1, . . . , 20 (see Figure EJ). 



FIGURE 7. Graph showing log 10 A^ 2N ' N \lj*) against N for N = 1, 
Thus for iV = 20, there are unit L^norm functions spanning the W. 

,(2N,N) 



,20. 



(2N,N) 



(era;*) 

and 7 3 (w*)>V2' :JV ' JV ' , ( | j-'*) subspaces whose difference in L 1 -nonn is less than 10~ 8 . 

Recall that for the cocycle in Theorem 2, the Oseledets subspace W2(w) in fact 
independent of uJi for i ^ 0. This contrasts with the current example: to see that 
here the Oseledets spaces W2(cj) do not depend only on ojq, it is enough to observe, 
for example, that ujq — uji = 1 but w^ N ' N \^*) and w^ N ' N \ct 2 lu*) are markedly 
dissimilar. 



8. Proof of the Multiplicative Ergodic Theorem for non-invertible 

matrices 

In this section we present a strengthened version of the Multiplicative Ergodic 
Theorem (MET) for the case of non-invertible matrices. Let a be an invertible 
measure-preserving transformation of (Q,B,^i) and consider a linear cocycle P : 
Z + x ft x K d — * M. d . Note that even though the matrices may be non-invertible, the 
invertibility of a is crucial to the argument. If the matrices are invertible then the 
two-sided cocycle is naturally defined as a map P : Z x ft x R d — ► M. d . 

Recall that in the case of a one-sided linear cocycle (satisfying certain mtegrabil- 
ity conditions) , the MET provides an invariant flag of subspaces of K d characterising 
the exponential growth rates of all vectors. For a two-sided cocycle, one obtains 
an invariant splitting of R d into Oseledets spaces by considering the intersection of 
each subspace in the flag of the forward cocycle with the corresponding subspace 
of the flag of the backward cocycle. Non-zero vectors v in the jth Oseledets space 
Wj((jj) satisfy Um n _ > ± 00 (l/n) log||P(ra, w,v)\\ -> Xj. 

In the case of a one-sided cocycle it clearly makes no sense to consider the limit 
lim n _ > _ 0o (l/n) log \\P(n, u>, v)\\ . Nevertheless one may still hope for an invariant 
splitting of R d rather than an invariant flag. This distinction is important if one is 
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interested in the vector corresponding to the one of the top characteristic exponents: 
the flag would only provide an invariant family of high- dimensional subspaces with 
the property that most vectors in the space have the correct expansion rate, whereas 
a splitting would provide an invariant family of low- dimensional subspaces, whose 
vectors are responsible for all expansion at the chosen rate. 

In this section we obtain a decomposition into Oseledets subspaces for a one- 
sided forward cocycle over an invertible transformation by means of a push-forward 
limit argument. 

Let || • || denote the matrix operator norm with respect to the Euclidean norm 
on R d . 

Theorem 3. Let a be an invertible ergodic measure-preserving transformation of 
the space (CI, B, /u). Let A: CI — ► M<j(R) be a measurable family of matrices satisfying 

J Iog + ||A(w)||d/i(w)<oo. 

Then there exist Xi > X 2 > ■ ■ ■ > Xt > — oo and dimensions mi,...,m<, with 

mi H h me = d, and a measurable family of subspaces Wj(u) C R d such that for 

almost every u> € CI the following hold: 

(1) dim Wj (u>) = mj ; 

(2) R d = ®U W ^)' 

(3) A(u>)Wj(u)) C Wj(auj) (with equality if Xj > -oo); 

(4) for all v £ Wj(w) \ {0}, one has 

lim (1/rj,) log || A(<r"— . . . A(c.;)t;|| = A,-. 

n — *oo 

Lemma 3. Let B: CI — ► M<j(R) be a measurable mapping into the space of symmet- 
ric matrices such that for almost all uo, B(u) has eigenvalues fj,i, . . . , m with multi- 
plicities mi,...,m<. Then there exists a measurable family ( e 3 Auj) ) 

V / l<j<t, l<i<mj 

of vectors such that the ^ej(w)^ form an orthonormal basis ofR d and e{(uj) lies in 
the iij eigenspace of B{uo). 

Proof. Consider the map R that takes a matrix and applies a single step of a row- 
reduction algorithm (e.g. find the first column that is not in row-reduced echelon 
form; transpose rows to put a non-zero entry in the correct place; divide so the 
leading coefficient is 1; subtract multiples of that row from all of the others; repeat) 
or does nothing in the case that the matrix is already in row-reduced echelon form. 
The domains of the pieces are measurable and therefore R is measurable. For all 
matrices A, R n {A) is a convergent sequence so the limit RRE(A) is a measurable 
function of the matrix. 

A collection of vectors spanning the kernel of a row-reduced matrix may be 
obtained in a measurable way. These vectors may then be measurably converted to 
an orthonormal set by applying the Gram-Schmidt orthogonalization algorithm. 

We apply this by taking a symmetric matrix B with eigenvalues with 
multiplicities mi,...,m<. We find an orthogonal set of vectors with each of the 
eigenvalues by applying the above procedures to B — /j,jL. Since all operations are 
measurable the proof is complete. □ 

Lemma 4. Let a: fl — > CI be an invertible ergodic measure-preserving transforma- 
tion and let (f n )^ =1 be a subadditive sequence of functions (that is a sequence such 
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that for every w G and each m and n, f n + m (uj) < f n (u>) + f m (a n uj)). Assume 
further that max(/i,0) is an L 1 function. Then there is a C G [—00,00) such that 
for almost every uj one has f n (ui)jn — > C and f n {o~~ n uj) jn — > C. 

Proof. The fact that there is a C such that f n /n — > C is Kingman's subadditive 
ergodic theorem. Letting g n (uj) = f n (a~ n u>), we see that g n + m (uj) < g n (u>) + 
gm,(o-~ n uj) so that the subadditive ergodic theorem applies to g n also (with the 
measure-preserving transformation being cr _1 ) and there is a constant D such that 
g n (uj)/n — > -D for almost all w. 

Since / n /n converges pointwise to C it also converges to C in measure. Similarly 
g n /n converges in measure to D. Since f n /n and g n /n have the same distribution, 
the constants to which they converge in measure must be equal. □ 

Lemma 5. Let a: Q ~* be an invertible ergodic measure-preserving transforma- 
tion and let A: Q. —> M<j(R) be a measurable family of matrices satisfying 

J log + \\A(u)\\ d/i(w) < 00. 

Let S be the multiset of characteristic exponents. Given uj S Q, let SV^(u>) be the 
multiset of logarithms of the nth roots of the singular values of A^ n \a~ n uj). Then 
for almost every uj, SV^(uj) -> S. 

Proof. Consider the family uj 1— ► A t (uj) with respect to the dynamical system 
er" 1 . Let the characteristic exponents be the multiset S'. This means that letting 
Sy'( n ) ( w ) be the multiset of nth roots of singular values of A T (a~~ n uj) . . . A T (o-~ x uj), 
one has SV'^ (u>) — > <S" for almost every ui S il by the MET. Since singular values 
are preserved by taking transposes we see that SV'^(uj) = SV^- 71 ' (uj). Thus it 
suffices to prove S — S'. To see this, note that (1/n) log || /\ k A^(uj)\\ converges to 
the sum of the largest k members of S, and (1/n) log || /\ k A 1 -(a~ n uS) . . . ^4 T (<7 -1 w)|| 
converges to the sum of the largest k members of 5", but these limits are equal by 
Lemma |4] □ 

Proof of Theorem [3J In the course of the proof we shall repeatedly use the symbol 
C to denote various constants depending only on uj. 

We write for the matrix product A(a n ^ 1 uj) . . . A(uj). From standard 

proofs of the MET, we have that [A^ (uj) t A^ n ' (uj)} 1 Z^ 2 ™' is convergent to a positive 
semi-definite matrix B(uj), for almost all uj, with eigenvalues e Al > . . . > e Xt with 
the correct multiplicities. We therefore let (el(uj)) be as in Lemma[3]and let Uj(ui) 
be the subspace of W 1 spanned by {e^(w): 1 < i < rrij}. The standard proofs of 
the MET show that if one lets Vj(ui) — ®\ = j Ui(ui) then the vector spaces Vj(u>) 
satisfy 

(1) A(w)Vj(u) C V^aco); 

(2) For all v € V s (u) \ V j+1 (cj), lim JWOC (l/n) log \\A^(uj)v\\ -> A,; 

For j < £, let W ( / l) (uj) = A^ n} (a~ n ui)Uj (o~~ n uj) and let W e (u) = U e (u). Then 
we claim the following: 

(1) For j < £, W^ n \uj) converges to an mj -dimensional subspace Wj(u>); 

(2) A(uj)Wj(uj) C Wj(auj); 

(3) If x S Wj(u) \ {0}, then (l/n)log||AW(o;)x|| -> Xj. 

(4) V j+1 («>)®W j (u>) = V j (u>). 
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Notice that W^ n+1 \auj) = A(u)W^ n \u) so that in the case j < £, ^ follows from 
(HJ) and the definition. For j = £, © and © follow from the standard MET proofs. 

Fix &j<£ and consider a basis Bq(lj) = {e|(w) : k > j, i < m^} for Vj+\(w) and 
a basis B\{lo) = {e* (w): i < rrij} for Uj(u>). The union of Bq(uj) and -Bi(u;) gives 
an orthonormal basis for Vj(ui). Since A(w)V 7+ i(w) C V, + i(ctcj) and A(uj)Vj(u>) C 
Vj (crw), it follows that if we express the linear transformation represented by A(w) 
with respect to the bases -Bi(w) U Bq(u>) and Bi(auj) U So(cra;), the matrix is of the 
form 



L ^ V^io(^) Ara (w) 

where if Vj+i(o;) is of dimension q = TOj+i + . . . + m^, the matrices Au(u>), Aio(cj) 
and Aqq(uj) have dimensions m_,- x rrij, q x rrij and q x q respectively. 

By definition, L^ n \u)) — L(a n ^ 1 uj) . . . L{u>). By analogy with the above we name 
the components of this matrix as follows: 



A$(aO A^( W ). 



We will need the following matrix identities: 


Claim 1. 


With A^ 


defined as above we have 


(8) 


a£V) 


= A n (a n - 1 uj)...A 11 (oj) 


(9) 




= A 00 ((J n ~ 1 uj)...A o(uj) 


(10) 




n-l 

= J2A$(a n - k u)A w (a 

k=0 



^)An '(w). 

Proof. The first two equalities are immediate and the third follows by induction on 
n. □ 

Claim 2. For almost every oj g J7, log ||Aoo^(w)|| — » Aj+i as n — > oo. 

Proof. One has for each i > j and 1 < fe < m !; (1/n) log |] A^ n \uj)e l k \\ — ► Aj, 
by the MET. It follows that considering A(™)(o>) as a linear map on Vj+i(uj), 
(l/n)log||A(«)(a;)| F . +l(a;) || - A, +1 . Thus 

(l/n)log||A^H||^A J+1 . 

□ 



Claim 3. .For every e > and for almost every w S fi, i/iere is D±(u>) such that 
1 oo 



||A^ l) (cr-"w)|| < J D 1 (w)e Tl ^+ 1+£ ) /or a// n > 0. 



Proof. Let f n {uj) = log ||Aqo (w)||. This is a sub-additive sequence of functions and 
f n (u)/n — > log Aj+i for almost every a; by Claim[51 Applying Lemma[2]we see that 
fn{o~~ n u)/n logAj+i for almost every ui. The claim follows. □ 

Claim 4. For every e > and for almost every oj £ Q, there is a F>2(w) < oo such 
that for alln>0 one has \\ A w {o~ n uj)\\ < D 2 (uj)e en . 
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Proof. By hypothesis log || is an integrable function and hence by a standard 

corollary of Birkhoff's theorem one has log || A(a~ n oj)\\/n — > 0. It follows that 
\\A(o-- n uj)\\ < D 2 (Lo)e en for a suitable D 2 (lo). Since ||Ao(w)|| < \\A(u)\\ the result 
follows. □ 



/ T \ 1/(2n) 

Claim 5. Under the above conditions, IA^\lu) A^l\u;)\ — ► e Xj I m . . 

Proof. To see this it is sufficient to show that every non-zero vector in Uj(u>) has 
growth rate Xj. Let u G Uj{uS) have expansion u — J2i< mj v i e i{ UJ )- 

First we show that A^(uj)v doesn't grow any faster than A^l' (u>)v. Note that 



\\A^(uj)u\\ 2 = i|A^ ) (cj) u || 2 + ||4o ) (^)w|| 2 sothatwehave||4o ) (^)«ll + ll4i ) (w)v|| 
grows at rate Xj . Applying the MET to A^ (u>), we see that A^ (u>)v grows at som 
rate A. We will show that A^q (u>)v grows at a rate no greater than max(A, Xj+i) 
It will follow that A = Xj . 
Equality (fTTJ)) gives 



\\A$(w)v\\ < E \\A^\a^ k w)\\\\A w (o n - k - 1 w)\\\\A^ k - 1 \u)v\\ 



k=0 

Fix an arbitrary e > 0. Claim |3] shows that \\A^ ((j n ~ k u;)\\ < D 1 (a n uj)e k ^+ 1+t \ 
Using Claim d] also, we see 

n-l 

l4oV>ll < D 1 {a n iu)D 2 {u J )Y / e k{X3+1+e) - ein ' k ' 1) - {A+e){n - k - 1) 



k=0 

There exists M such that -Di(w) < M on a positive measure subset of f2. By the 
ergodicity of a, there are infinitely many n for which D\ (er"w) < M. For these n, the 
right hand side of the inequality is bounded above by D 2 (w)Mne n ' mm ' Ai+1 ,A ' +2£ * . 
It follows that liminf log IIAjq^w)^!! 1 /™ < max(Ay+i , A) and thus 

limmflogll^tuOull 1 /™ < max(A j+ i, A). 

Since on the other hand limlog || A^ n \oj)u\\ l / n = Xj, we conclude that A > Xj. Since 
\\A^(lu)u\\ > we have Xj > A so that A = Xj as required. □ 



We now estimate 

q n (u!) = max — ^ , 

ves, \\ A W(a-»u)v\\ 

where Si denotes the unit sphere in M. m ' . Note that by scale-invariance one could 
equivalently define g n by taking the maximum over M. mj \ {0}. 
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We have 



g n (u>) — max 
veSi 





-^)ii ((T-( fc+i ^)4r^ 1) (^^)'' 









n-l 

< > max 

k=0 



n-l 

Emax 

fc=0 



(n-fc-1) , 



4i +1 V- (fc+l M? 



< 



n-l 

E 



max„ eS i 




min ueSl 


4i +1 V- (fc+i M« 





Note that in the third line we are making use of the fact that k ^ (a~ n uj) 
is invertible. 

Let e < (Aj — A 3 +i)/4 be fixed for the remainder of the proof. By Lemma[5]and 
Claim [5] the feth roots of the singular values of A^ (a~ k Lu) all converge to e Xj . It 
follows that there is a C > depending on to such that for all k, 



(11) 



mm 



> 



We remark that similar uniform lower bound estimates appear in the paper of 
Barreira and Silva [BS05j . Using Lemma [5l Claim [3] and Claim [4] there exists a C 
depending on ui such that for all k, 



max 



A^{a- k Lo)A w (a-( k+1 ^)i 



< C"e fc ^ Aj '+ 1+e )e efc . 



Combining the estimates we see 



n , n-l 

. 9n M<^E efc(Aj+1 ^ +3£) - 

fc=0 

Since 3e < Xj — Xj+t it follows that defining M(oj) = sup„ g n {w), one has M(oj) < oo 
for almost all u). 

We define a distance D between two subspaces of K d of the same dimension by 
the Hausdorff distance of their intersections with the unit ball B\ in M. d . We now 
estimate D (w^ n \tj), Wj m) (u;)) for m > n. 

Let x belong to the unit sphere of Wj (u>) (the distance is always maximized 
by points on the boundary). Then x = A^ {<7~ n uj)u for some u £ Uj(a~ n u;). 
Since for almost all u>, the matrix A^ n \o~ m uj) is invertible, there exists al- 
most surely a u' £ Uj{o~ m Lu) such that A^ m ~ n > (p~ m uS)v! — u + z where z £ 
Vj + i(o~ n Lu). Let v' be the coordinates of u' with respect to the basis Bi{uj). 



Then \\A { ™- n \a- m uj)v'\\ = \\z\\ and \\A { ™- n) (a 
\\z\\ < M{a- n Lj)\\u\\. Let y = A (m *> {(j- m Lj)u' so that y £ w\ m \uj). We then have 



LO )V 



It follows that 
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y = x + A( n \o~ n Lu)z. By Claim[3]we have 

\\A^ n \cj- n uj)z\\ < Ce ( ^+ 1+e) "||z|| 

(12) < Ce {x ^ 1+t)n M{<j- n Lo)\\u\\ 

for a C depending only on to. On the other hand, (llip implies that 

(13) 1 = \\x\\ = \\A^(a- n cu)v'\\ > C"e^- e ) n ||u|| 

for another C depending just on u). Let K = C/C and a = Xj — — 2e > 0. 
Dividing HU by (JTHJ) we see 



||y -a;|| = ||A(")(a- Tl w)z|| < Ke- an M{a-' n uo). 

The closest point of Wj m \u)) f] B\ to x is just the orthogonal projection of x onto 
Wj (ui) (which lies in Bi) so that the distance from x to Wj (u))nBx is bounded 
above by \\y — x\\ which in turn is bounded above by Ke~ an M(a~ n Lu). 

Conversely let y E B x n W^ m \u)). Then y = (cj- m uj)u' for some u' e 
Uj(a~ m ui). Let A'- m ~ n \<T~ m u))u' be decomposed into u + z with u <E Uj(a~ n ui) 
and z G Vj + \(a~ n uj). Let x = A( n )(<7 _n u>)it. Since sup n g„(w) = M(w), we have 
Ml < M(a- n u))\\u\\. So \\A {n \(j- n Lj)z\\ < KM(a- n uj)e~ an \\A { - n \(T~ n uj)u\\. We 
also have 

\\A^ n \a- n uj)u\\ < \\A^(^ n oj)(u + z)\\ + \\A (n \a- n Lu)z\\ 
(14) < l + KM(a- n ui)e- an \\A ( - n \a- n u!)u\\. 

So ||A( n )(o--"w)M|| < 1/(1- KM{a- n uj)e- an ), provided KM(<j- n u>)e~ an < 1. 
Combining this estimate with (fl4|) gives 



1 - KM(a~ n ui)e- 



As before it follows that the closest point of Wj (ui) n £?i to y is at a dis- 
tance at most KM(a~ n u))e~ an / '(1 — KM(a~ n io)e~ an ). In particular, provided 
that ifJV/((T- n w)e" Qn < 1, we have 



J ftTM(CT-"w)e- Q 
- ATM (a~ n ui)e~ 



fl(<H,^»H)< i 

Obviously for m,m' > n one then has 



£>( W" , (w),Ty (m ' ) (w) N ) < - 

Since M(uj) is measurable and tr is ergodic, there exist for almost all lo arbitrarily 
large values of n such that M(<j~ n us) < A for some fixed A > 0. It follows that the 
sequence of subspaces is Cauchy and hence convergent to a subspace Wj(u>). 

Let x belong to the unit sphere of Wj(u>). Then x — A^ n \a^ n Lo)u. As before, 
writing x as y + z with y G Uj(u>) and z € Vj+i(u>), we have ||z|| < M(w)||y|j. Since 
||j/|| 2 + ||z|| 2 = Xj we have ||y|| 2 (l + Af(w) 2 ) > 1 so that ||y|| > + M(u>) 2 =B. 

Thus each point of the unit sphere of Wj(u>) has a component in Uj(u>) of norm at 
least B. It follows that Vj(ui) = Vj+i(u>) © W^-(w). 

This completes the proof. □ 
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